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Abstract 

We study a class of nonlinear nonlocal cochlear models of the transmission line type, 
describing the motion of basilar membrane (BM) in the cochlea. They are damped dis- 
persive partial differential equations (PDEs) driven by time dependent boundary forcing 
due to the input sounds. The global well-posedness in time follows from energy esti- 
mates. Uniform bounds of solutions hold in case of bounded nonlinear damping. When 
the input sounds are multi- frequency tones, and the nonlinearity in the PDEs is cubic, we 
construct smooth quasi-periodic solutions (multi-tone solutions) in the weakly nonlinear 
regime, where new frequencies are generated due to nonlinear interaction. When the input 
is two tones at frequencies f\, (f\ < /2), and high enough intensities, numerical results 
illustrate the formation of combination tones at 2/i — /2 and 2/2 — fi, in agreement with 
hearing experiments. We visualize the frequency content of solutions through the FFT 
power spectral density of displacement at selected spatial locations on BM. 
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1 Introduction 



Digital signal processing on sounds is an essential component of modern hearing devices ^Ej , and 
a useful tool for evaluating acoustic theories of peripheral auditory systems, jTSI among others. 
A fundamental issue is to model the auditory response to complex tones because the nonlinear 
interaction of acoustic waves of different frequencies allows for audio compression ^B] among 
other applications. Nonlinearities are known to originate in the cochlea and are further modified 
in higher level auditory pathways. The cochlear mechanics has first principle descriptions, and 
so partial differential equations (PDEs) become a natural mathematical framework to initiate 
computation. However, in vivo cochlear dynamics is not a pure mechanical problem, and neural 
couplings are present to modify responses. To incorporate both aspects, a first-principle based 
PDE model was studied in [2T] for voice signal processing, where the neural aspect is introduced 
in the model phenomenologically. The first principle based PDE approach is more systematic 
compared with filter bank method ^H|, and has shown encouraging results. In |2*T] . time domain 
computation on multi-tone inputs revealed tonal suppressions in qualitative agreement with 
earlier neural experimental findings. 

In this paper, we shall analyze the well-posedness and construct multitone solutions of such 
PDE models in the form: 

p xx - Nu tt = e s (x)u t , xe(0,L), (1.1) 
p = mUft + r(x, \u\, \u t \)u t + s(x)u, (1.2) 

where p is the fluid pressure difference across the basilar membrane (BM), u the BM displace- 
ment, L the longitudinal length of BM; N a constant depending on fluid density and cochlear 
channel size; e a {x) > is the damping of longitudinal fluid motion; m, r, s are the mass, 
damping, and stiffness of BM per unit area, with m a constant, s a continuously differentiable 
nonnegative function of x. The coefficient r is a nonlinear function(al) of x, u, uf. 

r{x,\u\,\u t \) = r a (| Mi | 2 )+ 7 P{\u{x', t)\)K{x - x') dx' . (1.3) 

Jo 

Here: (HI) r a (-) is the local part of BM damping, it is a nonnegative continuously differentiable 
monotone increasing function, r a (0) > 0. In the nonlocal BM damping: (H2) K = K(x) is a 
localized Lipschitz continuous kernel function with total integral over i6i?' equal to 1; (H3) 
P(-) is a nonnegative continuously differentiable function such that for some constant C > 0: 

P(0) =0,P(q) < C(l + q 2 ), Wq> 0. (1.4) 
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The boundary and initial conditions of the system are: 

p x (0,t) = T MPT (t) = /(*), p(L,t) = 0, (1.5) 
u(x, 0) = uo(x), Ut(x, 0) = ui(x), (1.6) 

where the initial data is such that (u ,ui)(L) = (0,0); prif) is the input sound pressure at the 
eardrum; and Tm is a bounded linear map modeling functions of middle ear, with output de- 
pending on the frequency content of pr{t)- If Pt = J2j=i Aj exp{iujt} + c.c, a multi-tone input, 
c.c denoting complex conjugate, Jm a positive integer, then TmPt^) = Z)/=i Bj explicit} + c.c, 
where Bj = a M (ujj)Aj, c.c for complex conjugate, a^O) a scaling function built from the filtering 
characteristics of the middle ear [5]. 

Cochlear modeling has had a long history, and various linear models have been studied at 
length by analytical and numerical methods, jH], JO] and references therein. A brief derivation 
of the cochlear model of the transmission line type, e.g. the linear portion of (jl.ljl - fjl.fijl . is 
nicely presented in [TH] based on fluid and elasticity equations. 

It has been realized that nonlinearity is essential for multitone interactions, [SJ 13 12 H] etc. 
Nonlinearity could be introduced phenomenologically based on spreading of electrical and neural 
activities between hair cells at different BM locations suggested by experimental data, jjj, j^j. 
Such a treatment turned out to be efficient for signal processing purpose [2T], and (jl.3j) - (jl.4|) is 
a generalization of existing nonlinearities [7], j^j, [TO] . 

Multitone solutions require one to perform numerical computation in the time domain. The 
model system (jl.l)) - (jl.6j) is dispersive, and long waves tend to propagate with little decay from 
entrance point x = (stapes) to the exit x = L (helicotrama). The function e s (x) is supported 
near x = L, its role in numerics is to suck out the long waves accumulating near the exit |21j . 
Selective positive or negative damping has been a novel way to filter images in PDE method of 
image processing jT^j. In analysis of model solutions that concern mainly with interior properties 
however, we shall set e s to zero for technical convenience. 

The rest of the paper is organized as follows. In section 2, we perform energy estimates of 
solutions for the model system (jl.lj) - ()1.6|l . prove the global well-posedness and obtain growth 
and uniform bounds in Sobolev spaces. In section 3, we construct exact multi-frequency solu- 
tions when 7 is small enough and nonlinearity is cubic, using contraction mapping in a suitable 
Banach space. The constructed solutions contain all linear integral combinations of input fre- 
quencies. In section 4, for two input tones with frequencies fi and f 2 (fi < f 2 ), we illustrate 
numerically the generated combination tones 2fi — f 2 and 2f 2 — j\ on power spectral density 
plots at selected points on BM. These tones are heard on musical instruments (piano and violin), 
in particular, 2fi — f 2 is known as the Tartini tone. The conclusions are in section 5. 
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2 Global Well-Posedness and Estimates 

Let us consider the initial boundary value problem (IBVP) formed by and show that 

solutions exist uniquely in a proper function space for all time. To this end, it is convenient to 
work with the equivalent integral form of the equations. It follows from (jl.lj) and (jl.5|) that: 



.I- 



Vx= \ (Nu tt + e s (x)u t )dx + f(t), 



o 



L 



-p{t,x)= dx' {Nuu + e s u t )dx" + f(t)(L-x). (2.1) 

Jx JO 

Combining (JOJ) and (|Ol . we get: 

mu tt + [ dx' [ (Nu tt + e s u t ) dx" + f{t)(L - x) = -r(x,\u\,\u t \)u t - s(x)u, (2.2) 



o 



with initial data (jl.6|) . Let w = (wi,w 2 ) = {u,u t ), and write ()2.2|) into the system form: 

w ltt = w 2 , (2.3) 

pL px' pL px' 

mw2 1 + / dx' / iV W2,t dx" = — dx' / e s w 2 dx" 
Jx Jo Jx Jo 

—r(x, \wi\, \w 2 \) w 2 - s(x) wi + f(t) (x - L). (2.4) 

The related integral form is: 

Wi = u + [ w 2 (x,r)dr, 

pt pL 



o 



Pt PL, PX 

Aw 2 = Au\ — dr dx' / e s w 2 dx" 
Jo Jx Jo 



t pt pt 

drr(x, \wi |, \w 2 \) w 2 — s(x) I widr + (x — L) f(r)dr, (2.5) 
o Jo Jo 

where A : L 2 ([0,L]) — > L 2 ([0,L}) is a bounded self-adjoint linear operator: 

rL px 1 _, 

Ag = mg + dx' N g dx" = m g + A g. (2.6) 

Jx Jo 

To see the self-adjointness of A, let g, h G £ 2 ([0, L}), then (Ag, h) = (g, Ah) or: 

dx' [ X g{x")dx") h{x)dx= [ L \ t dx' T g(x")dx") d [* h 

\Jx Jo ) Jo \Jx Jo / JO 



o \Jo 



dx' g(x')^j (^J dx' h(x') \ dx, 
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hence (Ag, h) L 2 = 
square L 2 norm: 



(g,Ah) L 2. Clearly, A is bounded; also (A-,-) = (A-,-) L 2 is an equivalent 



j-L rx 

(Ag,g) = m\\g\\ 2 2 + N ( gf dx > m\\g\\ 2 2 , 

Jo Jo 

(Ag, g) < m\\g\\l + NL 2 \\g\\ 2 2 = (m + NL 2 )\\g\ 



(2.7) 



Moreover, A has a bounded inverse. To see this, note that A — m Id is a compact operator on 
L 2 ([0, L]), so Riesz-Schauder theory |22] sa y s that the spectrum of A can have only eigenvalues 
of finite multiplicities except at number m. On the other hand, zero cannot be an eigenvalue of 
A, as (Ag,g) > rn\\g\\ 2 . The bounded inverse of A follows, and we denote it by A^ 1 below. 

Now we establish the global existence of solutions of ()2.3j) - (J2.4j) in the function space 
C([0, oo); (H l ([0, L])) 2 ). It is straightforward to show by contraction mapping principle that if 
\\(u , Ui)\\ H i < oo, there is a time £* such that ()2.5|) has a unique solution in C([0, £*); (i/ 1 ([0, £])) 2 ) 
under our assumptions on the nonlinearities. Such a solution in fact lies in C 1 ([0, £*); (if 1 ([0, -^])) 2 ), 
and obeys the differential form of equations (|2.3|) - (j2.4j) . with both sides interpreted in the H 1 
sense. Taking the limit x — > L, we find that the system ()2.3|) - ()2.4j) reduces to the ODE system: 



with initial data (wi, w 2 )(L, 0) = (0,0), hence (wi,w 2 )(L,t) = (0,0), Vt G (0,t*). 

Let us derive global in time estimates of solutions in H 1 to extend the local solutions to 
global ones (so t* = oo). The left hand side of (|2.4|) is just (Aw 2 )t, and: 



w 2 ,t = -r{t) w 2 - s(L)w!, 



(w 2 ,(Aw 2 ) t ) = (w 2 ,A(w 2 ) t ) = (Aw 2 ,w 2>t ) = (w 2)t ,Aw 2 ), 



so: 



d_ 

dt 



(Aw 2 ,w 2 ) = (w 2>t ,Aw 2 ) + (Aw 2 j,w 2 ) = 2(w 2 ,Aw 2 ,t), 



(2.8) 



hence ~^(Aw 2 ,w 2 ) = (w 2 , Aw 2)t ). Multiplying ()2.3|) by Wi, f!2.4|) by w 2 , adding the two expres- 
sions and integrating over [0, L], we estimate with Cauchy-Schwarz inequality: 




+ (w 2 , wt) - (sw u w 2 ) + (fit) (x - L),w 2 ) 



< -r a (0)|| W2 || 2 + r a (0)||«; 2 || 2 + -J—\\f(t)(x - L)|| 2 
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^ (j X e s w^j (f^) dar+||(l-s)||oo|k2||2|K||2 

< \f(t)\ 2 L 3 + L a ||e.|| 2 |MI2 + - s\U\\ Wl \\l + \\w 2 \\l) 

< (~l|i - «IU + £ 2 || es || 2 )|K|| 2 + i||i - -HoolKlll + |/(t)| 2 ^ 3 . (2.9) 

Let Ci = max(^||l-s|| 0O + 2 ^||e s || 2 ,|||l-s|| 0O ), and 2E = (w^Wi) + (Aw 2 ,w 2 ), we have from 
(EH): 



or: 



E(t)<E(0) + C 1 ( E{s)ds + —— [ \f\\s')ds'. 
Jo 12r a {0) Jo 

Gronwall inequality implies: 



£(i)£ ( £(0) + l2S)/o' l/|2(iVi ') eC '' 



or: 



||K,^ 2 )|| 2 <min(l,m)- 1 ^(0) + T ^^/ o ' \f\ 2 {t')dt'^ e c « . (2.11) 

Next we obtain the gradient estimates. Differentiating ()2.3|) - (j2.4j) in x gives: 
d 

—w 1>x = w 2>x (2.12) 

^- (mw 2 x — N ( w 2 (x', •) dx'^] = ( eJx') w 2 {x') dx' — rw 2 x — s'w 1 — swi x + fit) 
at \ Jo J Jo 

- \2r a w 2 w 2tX + f L P(w 1 )(x',t)K x {x-x')dx' \ w 2 . (2.13) 



o 



Multiplying (|2.12|) and (|2.13|) by Wi jX and w 2 ^ x , and integrating over x G [0, L], we find: 
1 d ,,, ,, n „ ti /" L 



Id d f L f x 

o3t(II w i,x|| 2 + ^lk 2ja; || 2 ) = iV— / w 2jX / iu 2 (a:', •) cfx' + iw 2i!B ) 
2 dt ~ dt Jo Jo 

^2r' a w 2 w 2>x + P(w x )K x \ w 2 w 2iX 



{u>2,x,J o e s w 2 ^j - rw\ x 



- s'wxW 2iX - / sw 1:X w 2:X + f(t) w %x . (2.14) 
Jo Jo Jo 

The integral in the first term of the right hand side of (|2.14|) equals: 

— dxw 2tX (x,-) / w 2 (x', •) dx' = -— / wl(x',t)dx', (2.15) 
2 Jo Jo 2 Jo 







where we applied integration by parts once and W2(L,t) = 0. 
The other terms are estimated as follows: 

dx P(wi)(x' 1 •) K x (x — x') dx'^j w 2 w 2jX 

< C f L \w 2 w 2iX \ [ L (l + \ Wl \ 2 )\K x \(x-x')dx' 

Jo Jo 

< C{l+\\ Wl g) [ L \w 2 w 2>x \ 

Jo 

< C(l+\\ Wl \\l)\\w 2 \\ 2 \\w 2 , x \\ 2 

< So\\w 2 Al + ^C\l + \\w4l) 2 \\w 2 \\l, (2.16) 

4d 

for any 5 > 0, C = Cdl-K^Hoo), by (jl.4|) . Integration by parts and w 2 (L,t) = give: 

px pL 

(w2,x, e s w 2 ) = - e s wldx<0. (2.17) 
Jo Jo 

Estimate with Cauchy-Schwarz inequalities to get: 

L rw\ x - 2 [ L r» 2 % < -r a (0)||w 2)X ||l; (2.18) 
Jo 

- s' w x w 2x dx- s w lx w 2tX + (w lx , w 2x ) 
Jo Jo 

< ||1 - sHoollwi^l^lluvlh + ||s'||oo||Wi||2||w 2 ^||2 

< ^o\\w 2 ,x\\l + y^HIi - sWU^iAl + ^HKirJKili; (2.19) 

/(f) [ L w 2 , x < \f{t)\L^\\w 2 , x \\ 2 < 5 \\w 2 , x \\ 2 2 + -^\f(t)\ 2 L. (2.20) 
Jo 4oo 

Combining (j2~T4ll - (l2~2711) . with A5 = r a (0)/2, we get: 

^(IK^I* + HI«**II2) < -~ IMI2 - ^II^Hl + C 2 (t) + C 3 (t)\\ Wl , x \\l (2.21) 
where: 

C 2 (t) = + MM^\\ Wl \\l+ ^(1 + \\ Wl \\ 2 2 ) 2 \\w 2 \\l, (2.22) 
4o looo 4o 

C3(*) = ^ r l|l-*. (2-23) 
16d 

and ||(u>i, w 2 \\ 2 are bounded as in (J^TTTJ). Integrating ()2.2H) over t G [0,T], we find: 

N r T r T 

m , \\(w h x,w 2 ,x)\\l(T) + -\\w 2 \\l<C A + C 2 (t')dt'+ C 3 (t')\\(wi,w 2 ) x \\l, (2.24) 

2 Jo Jo 
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where m! = | min(l, m), 

C 4 = ^max(l,rn)||(w , :i; ,wi, a; )||2 + y 1 1 111; (2.25) 
and Gronwall inequality implies: 

HK„^)||2(T)<— C 4 + / C 2 (t)dt exp{— / C 3 (t)dt}. (2.26) 
T7i y jo j m Jo 

We see from (J2H that io 2 , t G C{[0,oo); H\[0, L})), hence pressure p G C([0, 00); F 3 ([0, L])) 
from (|2.ip . We have thus shown: 

Theorem 2.1 Under the growth condition \1.4\) o,nd the initial boundary conditions il., 5j) and 
< li.6)) . i/ie model cochlear system 4i.i|) -i f?~3)) /ias unique global solutions: 

(u,u t ,p) G C([0,oo); (^([0,L])) 2 x # 3 ([0,L])). 

The estimates can be improved with the additional assumption: 

3r (0) 

s{x)>s o >0,Vxe[0,L\; ||ej 2 < — (2.27) 

2L2 



Theorem 2.2 (Growth Bounds) Under the additional assumption \2.21 ), the global solu- 
tions in Theorem 2.1 satisfy the bounds: 



(u,u t )\\l+ f \\u t \\l(t>)dt> < ai + a 2 f \f(f)\ 2 dt', 
Jo Jo 

i-t rt rt' 



11(^,^)111+/ |Kt||!(0 dt' < a 3 + aj (1+/ \f(t")\ 2 dt"fdt\ (2.28) 
Jo Jo Jo 

for some positive constants ai, i = 1,2, 3, 4. 

Proof: Multiplying (|2.3j) by s(x) Wi, and (|2.4j) by w 2 , adding the two expressions and integrating 
over [0,L], we estimate with Cauchy-Schwarz inequality: 

(s Wx, w ht ) + (w 2 , (Aw 2 ) t ) = - i^J dx ' J Q e sW 2 dx", w 2 ^ 
-(r w 2 , w 2 ) + (f(t) (x-L), w 2 ) 
< -( w 2 (x') dx', f e s {x')w 2 {x') dx 1 ) - r a (0)\\w 2 \\ 2 2 + f(t)(w 2 , (x - L)) 



1 I f I 2 

< -(r.(0) - -L^\\e s \\ 2 )\\w 2 \\\ + 5\\w 2 \\l + ^\\{x - L)\\\. (2.29) 
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Choose 25 = r a (0) - |L 3 / 2 ||e s || 2 > to find: 

Id I f I 2 L 3 

+ (Aw 2 ,w 2 )) < -5\\w 2 \\l + ^-y. (2.30) 

So: 

J((swiK + (A™ 2 ,™ 2 ))(t) < -<5 Alalia +co + ci f \f(t')\ 2 dt'. (2.31) 
2 jo Jo 

where c = |((sti>i, Wi) + (Aw 2 ,w 2 ))(0), and ci = Hence, 

min(sn, m) .. , , ,,o f',, ,,„ ,, 

-iik^ < c + Cl / i/(or< 



Z JO 



/ ||u; 2 ||'(t')^' < r^o + r 1 / |m')l 2 ^'- (2-32) 

JO JO 

In particular, if /(£) a bounded continuous function, (|2.32j) gives the growth bounds: 

||(™i,OI| 2 < 0(^2), /' || W2 || 2 (t')^<0(t), (2.33) 



o 

implying that ||u> 2 || 2 has a bounded time averaged L 2 norm square. 

Similarly, we improve the gradient estimates. Multiplying ()2.12|) by swi^ x , (j2.13J) by w 2tX , 
integrating over x G [0, L], we cancel out the two integrals on s(x) Wi <x w 2<x . Proceeding as 
before, we arrive at: 

dl 2 N d ll2 r a (0) .. Il2 , N 
Tgll^M^M) + r "||w 2 , a; || 2 ) < ~Ydi^ W2 ^ 2 iFa^-L + C 2 (t), (2.34) 

and so integrating over [0, T] gives: 

-((swi tX ,w 1>x ) + m\\w 2 , x \\l)(T) + / H^^ll^t')^' + irlk 2 || 2 (T) 

z z Jo 2. 

1 N f T f T 

< -((su , x ,u , x ) +m\\u hx \\l) + — |K||1+ / C 2 {t') dt' = C 5fi + C 2 {t')dt'. 

2 Z JO JO 

(2.35) 

If s(x) > s > 0, then: 

UK*, «> 2> s) lll(r) + ^ ||^||l(t) dt < ^ fc 5i0 + jf C 2 (t) dtj , (2.36) 

where m" = | min(s , m). Substituting ()2.32|) in ()2.36|) gives (|2.28|) . In particular, for a bounded 
continuous /(£), 

HK^^II^T) + jf H^IIKO tft < 0(T 3 ). (2.37) 
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The proof is finished. 

If the nonlinear damping functions are bounded, i.e: 



r„(0 < C 6 , < C 6 , P(0 < C 6) V£ > 0, (2.38) 

for some positive constant C*6, then we have: 



Theorem 2.3 (Uniform Bounds) Under the assumptions \2.21\) and \2.3$i) . and that f(t) 



is a bounded continuous function, the global solutions in Theorem 2.1 are uniformly bounded: 

\\(u,u t )\\ H i + \\p\\h* < C 7 < oo, Vt > 0, 
for some positive constant C\ . Moreover, the dynamics admit an absorbing ball: 

limsup(||(M,Wt)|| H i + IMIfl 3 ) < c s, 

t—*oo 

where Cg is independent of initial data. 

See JE] for an example of a bounded damping function. The energy inequality ()2.3())1 lacks a 
term like —const. H^iHl on the right hand side, and so is insufficient to provide uniform bounds. 
The idea is to bring out the skew symmetric part of the system. 

Proof: multiply ()2.3|) by mw 2 , ()2.4|) by W\, integrate over x G [0, L], and add the resulting 
expressions to get: 

m(w 1 ,w 2 ) t + (Aw 2 ,t, wi) = HIW2II2 _ ( sw u w i) - ( rw 2, 

+ (f(t)(x — L), Wi) - (w u [ L dx' f e s w 2 dx"). (2.39) 

Jx JO 

Using the identity: 

■^(Aw 2 ,w 1 ) = (Aw 2) t,w 1 ) + (Aw 2 ,w 1;t ) = (Aw 2)t ,wi) + (Aw 2 ,w 2 ), 

we have: 

d 



m(wi,w 2 ) + (Aw 2 ,Wi)] = (Aw 2 ,w 2 ) + m\\w 2 \\l - (sw l ,w 1 ) - (rw 2 ,wi) 



dt 

/•L rx 1 

+ (f(t)(x-L),w 1 )-(w 1 , dx' / e s w 2 dx") 

Jx JO 



< (NL 2 + m)\\w 2 \\l - s \\wi\\l + Ce 1 1 ^2 1 1 2 1 1 ^1 1 1 2 

+ \f\L 3 / 2 \\ Wl \\ 2 + L 3//2 ||e 2 || 2 |K||2 IK|| 2 . (2.40) 
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Apply Cauchy-Schwarz to polarize the last three terms to get: 

r< II II II II ^ I, ||2 . s n i|2 

so 4 
\f\L^\\ Wl \\2 < |/| 2 L 3 / So + ||K||^ 

L 3/2 \\e 2 \\ 2 \\wx\\ 2 \\w 2 \\ 2 < L3||e2|l ' ||w 2 ||^ + ^||w 1 || 2 , (2.41) 

it follows from ijCTTjl that: 
d 



L m(«;i,twa) + (Aw 2 ,wx)} < C 9 \\w 2 \\ 2 2 - jKlll + 1/1^7*0, (2.42) 



where: 



r 12 r 3 u ii 2 

-2 i i °6 , ^ ll e 2 || 2 



C 9 = iVL + m + — + 



so s 

Multiplying (|2.30|) by a positive constant C p and adding the resulting inequality to ([2.42)1 . 
we find: 

±E, < (-5C P + C 9 )\\w 2 \\l - ||M| 2 + |/| 2 + (2.43) 

where: 

C 

E p = -^((swx,wx) + (Aw 2 ,w 2 )) + m (w l7 w 2 ) + (Aw 2 ,wx). 
Choose C p large enough so that C p > Cg/S, and: 

C 

E p > -?-{{swx,Wx) + (Aw 2 ,w 2 )) 

> min( So ,m)^||( Wl ,«; 2 )|| 2 ; (2.44) 

so: 

3C 3C 
Ep < —^-{(swx,wx) + {Aw 2 ,w 2 )) < —2-max.(\\s\\ 0O ,m + iVL 2 )||(wi, w 2 )||2, 

thanks to (|2.7)h then for some positive constant Cxo- 

± Ep <-CxoE p + \f\\- + ^). (2.45) 
at s 12t> 

The uniform bound on || (u, u t ) \\ 2 follows from (|2.45j) . Moreover, the fact that Cxo is indepen- 
dent of initial data implies the absorbing ball property of (u, u t ) in L 2 (i.e. the limsup as t — > oo 
is bounded independent of initial data). Equation ()2.4|) and L 2 invertibility of the operator A 
imply a similar uniform bound on u tt - Equation (|2.1j) in turn shows that ||p||# 2 is uniformly 
bounded and has the absorbing ball property as well. 
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Now we proceed with gradient estimate of (wi,w 2 ). The symmetric inequality is just (|2.34j) 

but with C 2 uniform in time now. The skew symmetric inequality is obtained by multiplying 

mw 2) x to (|2.12j) plus Wi jX times ()2.13|) . integrating over x G [0, L\. 

d d f x 

— [m(w ltX , w 2 , x )} - N(w 1)X , — w 2 (x',t) dx) 



< (w lx , / e s w 2 ) - (rw 2x ,wi x ) - (s'wi,w lx ) - s \\w liX \\ 2 
Jo 

+f(t)(l, w ltX ) + 2C 6 (\w liX \, \w 2;X \) + C n (\w 2 \, \w 1;X \). (2.46) 



The second term on the left hand side equals: 



N(w 1 ,w 2>t ) = N^(wi,w 2 ) - iV||u> 2 ||2- 



It follows that: 



-j-[m(w 1>Xl w 2>x ) + N(w 1 ,w 2 )} < 4C 6 1| w x>x \\ 2 \\w 2>x \\ 2 - y \\w ljX \\l + C 12 , (2.47) 

for a positive constant C\ 2 ] or: 

d 16C 2 s 

— [m(w 1>x ,w 2tX ) + N(w 1} w 2 )] < -11^2,^111 — ^11^1,^111 + C12- (2.48) 

dt s 4 

Multiplying a constant C' p > to ()2.34|) with C 2 constant, and adding the resulting inequality 

to (l2~3Ki we get: 

jE p < -j\\wi,x\\ 2 + Cis - — )lF2,*|| 2 > ( 2 - 49 ) 

where: 



C 

E' p = -^((s^jKJi,!) -h 7V||-iJ7 2 1|| + mll^Ha) + m(w hx , w 2tX ) + N(wi, w 2 ), 

and C*i3 = C 2 C' p +Ci 2 . The term (wi,w 2 ) is bounded from above by constant times ||(iui,x, tf2,a;)||| 
due to Poincare inequality and so we can choose: 

32C| 



CL > 



-6 



P s r a (0Y 

and large enough so that for some positive constants Cu, C' l4 : 

C' u \\(w 1:X ,w 2>x )\\l > E' p > CuUw hx ,w 2 , x )\\ 2 2 . (2.50) 

Inequality (|2~4T)j) yields: 

j t E' p <-C 15 E' p + C 13 , (2.51) 

implying the uniform estimate on \\(u x , Ut tX )\\ 2 and the absorbing ball property. The uniform 
estimate and absorbing ball property on ||p||#3 follows from ()2.4j) and (|2.1j) . The proof is 
complete. 
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Remark 2.1 The estimates in Theorem 2.3 imply that the evolution map denoted by S(t) 
is relatively compact in the space (u,u t ) G (L 2 ([0, L])) 2 . Hence for any bounded initial data 
(u, u t )(0) G (i? 1 ([0, L})) 2 , the dynamics (u,u t ) approach, in the space (L 2 ([0, L])) 2 ; the univer- 
sal attractor A defined as: 

A = n t>Q S(t)B P0 , 

where B Po denotes the ball of radius po in (H l ([0, L])) 2 , the absorbing ball given by the estimates 
of Theorem 2.3. 



3 Multitone Solutions 

In this section, we consider special solutions to (jl.lj) - (jl.2|) that exhibit explicitly their frequency 
contents. For simplicity, let us assume that e 8 (x) = 0, and s(x) > s > 0, Vx G [0, L]. 

3.1 Linear Waves 

First we consider the linear regime with r = r , a positive constant. Solutions are superpositions 
of single frequency time harmonic waves of the form p = P(x)e lu)t + c.c, u = U(x)e luJt + c.c, c.c 
complex conjugate, where P and U are complex functions that satisfy: 

P = (-muj 2 + ir uj + s(x))U, (3.1) 

P xx + Nu 2 U = 0, (3.2) 

P x (0) = P in , P{L) = 0. (3.3) 

Let: 



—mu 2 + ir u + s(x) ' 
so: 

(a,(3) = Nuo 2 (-mu 2 + s(x), -r u)/[(-mu 2 + s(x)) 2 + r 2 Q u\ (3.4) 
then (j3.1|) - ()3.3|) is equivalent to: 

P xx + (a(x) + i/3{x)) P = 0, (3.5) 
subject to (USD. If uo + 0, (3 + 0. We show: 

Lemma 3.1 The boundary value problem, \3. 5)) and \3. 3]) . has a unique solution for alluo such 
that ||-P||_H-2([ 0ii ]) < C|Pj n |, for some constant C independent of to. 
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Proof: Write P = P- m {x — L) + Q, then Q satisfies: 

Qxx + (a(x) + i(3{x)) Q = ~{a{x) + i(3{x)){x - L)P in = {f x + if 2 ){x), (3.6) 

with boundary conditions Q x (0) = 0, Q(L) = 0. The left hand side is a Fredholm operator on 
Q, so it is sufficient to prove that zero is not an eigenvalue, which follows from the estimate 
below. Write Q = qi + iq 2 , then: 

(:;Hrf)(:;H£)- 

Multiplying ()3.7|) by (q 2 , —qi), we find: 

qi,xxQ2 ~ qiq2,xx ~ Pql - (3q\ = f x q 2 - f 2 qi, 
which gives upon integrating over x € [0, L], integrating by parts: 

- F P (l! + <&) dx = f\hq 2 - fm) dx. (3.8) 
Jo Jo 

Multiplying (|3.7|) by (qi,q 2 ), integrating over x G [0, L], we have after integration by parts: 

-11(^,^)112+ j L oc{ql + ql) = f L (hqi + f 2 q 2 ). (3.9) 

JO JO 

It follows from ()3.9j) and Poincare inequality that: 

||(?l,?2)||l < L 2 \\(q 1<x ,q 2 , x )\\l 

< L a ||a|| oe ||(gi l « & )||l + L a ||(/i,/ 2 )||2||(gi,82)||2. (3+0) 

For \u\ < oj < 1, Halloo < 2Nu 2 /s , 2lu 2 L 2 N/s < 1/2, 

||(gi,g 2 )|| 2 <2L 2 ||(/ 1 ,/ 2 )|| 2 . (3.11) 

As \\(h,f 2 )h = 0(u 2 ), (EUD implies ||(gi,g 2 )|| 2 = 0{u 2 ) for u « 1. For u 2 > M = 
M(7V, m, 5, r ), M large enough, a ~ — — , f|3.9jl shows: 

ll(?l,?2)|| 2 <^||(A,/ 2 )||2. (3.12) 
When uj 2 < uo 2 < M, |/3| is bounded from below uniformly in u: 



N r \u 



3 



1^1 — Ti 2~> TTU2 1 2~~2 — 

|| — mar + s^xJII^ + TqUJ z 
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for some positive constant (3q only depending on r , m, L, and s(x). Inequality (|3.8|) gives: 

\\(qi,q2)h < P 1 \\(fiJ2)h = Po'Cii^m,^^^)^ (3.13) 

for positive constant Ci(N, m, r , L, s), uniformly in u 2 G [uq, M]. Combining ()3.11|) . ()3.12|) . and 
(|3.13|) . we see that for any u, and any given P in , there is a unique solution P, \\P\\2 < C^l-fml, 
for constant Ci independent of lu. The lemma is proved by applying the L 2 estimate and the P 
equation (|3.5p . 

3.2 Nonlinear Waves 

We are interested in the persistence of multitone solutions when nonlinearities are present. For 
simplicity, we shall consider: (Al) r a > 0, a constant, and P(u) = u 2 , the overall nonlinearity 
is cubic. As for linear waves, assume that (A2) e s = 0, s(x) > s > 0, s G C 1 ([0, L]). We prove: 



Theorem 3.1 (Existence and Uniqueness of Multitone Solutions) Let the left bound- 
ary condition be: 

fin(t) = a i ex P{^} + c - c -, 

j=l,---,m 

and fix p > 1. Then under (Al-(A2) and for 7 small enough (independent of p), system M.l\) - 
has a unique solution of the form: 

u(x,t) = Uk(x) exp{ik ■ out} + c.c, (3-14) 

kez m 

where uo = (ui, 002, • • • , uj m ), and complex valued functions Uk{x) G H l ([0, L]), such that: 

\\u\\ =Y, P |fc| \\Uk\\m < 00. (3.15) 
k 

The pressure p is similar. 

Proof: Let B be the Banach space consisting of space-time functions of the form (|3.14j) with norm 
(|3.15j) . Let B>\ = {v G B : Vt G B}, a subspace of B. Consider the mapping M : v G B\ — > u 
defined as the unique bounded solution of the following equation in B\\ 

muu + / dx [ dx" Nu u + fin(t)(L - x) 
Jx Jo 

+r a u t + s(x)u = -jr n i(x,v 2 )v t , (3.16) 
where 7 r n i is the nonlinear nonlocal part of the damping function. 
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Let us show that M is a well-defined bounded mapping from B\ to itself. First we notice 
that for any functions Ui G B, i = 1,2: 



= E( E ^^,)e^', 

jeZ m k 1 +k 2 =j 

so: 

IKu 2 || = E II E ^iMfeall// 1 

jeZ m ki+k 2 =j 

< £ p^l + l^l|| MfclMfc2 || H1 

= E p'^'IKJk 1 E p'^WIIh 1 = IKHIMI- (3-i7) 

It follows that: 

IM^ 2 MII <c||v 2 *#IINI < cVIINI = c"lkH 2 IKII: ( 3 - 18 ) 

where C" depends on the kernel function K, and * denotes the convolution integral on [0, L\. 
Denoting F(x,t) = —^r n i(x,v 2 )v t G B, we show that u G B\. Write F(x,t) = Y^kez m Fk e lfc ' w * 
+ c.c, then (|3.16|) is same as the system: 

p xx - Nu tt = 0, 
p = mu tt + r a u t + s(x)u + F, (3.19) 

with boundary condition: p x (0, t) = fi n (t), p(L,t) = 0. Seek solution of system (|3.19J) in the 
form: 

P= E Pke^ + cc, 
kez m 

u = Wfce ifc ""* + c.c, 

k&Z m 

resulting in (k = [k\, k 2 , ■ ■ ■ , k m ))\ 

p k = {-m(k ■ to) 2 + ir a (k ■ u) + s(x))u k + F k , 

Pk, xx + N(k ■ ufu k = 0, (3.20) 

with boundary condition: Pk(L) = 0, Pk,x{0) = if k is not one of the m modes (along ej, 
j = 1, 2, • • • , m) of / in ; otherwise Pk,x{0) = if = e,-. 
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For each k, the system (13.20)) can be uniquely solved as in Lemma 4.1, with the estimates: 

\\Pk\\m < Ci||F fc || 2 + C 2 ||/i n , fc || 2 , 

IKIU* < C 3 \\ Pk \\ H i/(l + \k-w\). (3.21) 

It follows that the mapping M is from B\ to itself, and it is not hard to check that M is a 
contraction mapping if 7 is small enough. The proof is done. 

4 Numerical Results 

The model system is computed with a second order semi-implicit finite difference method, 
we refer to [21] for its details and choice of the coefficient functions in the model other than 
P(u) = u 2 . The input at the left boundary x = is the sum of two tones (sinusoids) at 
frequencies fx = 3.5 kHz (kilo Hertz) and / 2 = 4 kHz, with amplitudes 80 dB (decibel) and 
85 dB respectively. The zero decible is 20 \x Pa in physical unit. The time step is 0.01 ms 
(millisecond), and spatial grid is 0.01 cm. The computation ends at 20 ms when the BM 
responses reach a steadily oscillating state. To observe the frequency content of such a state, we 
select four points (x/s, j = 1, 2, 3, 4) on BM, and examine the response time series (u(xj, t)) at 
these points from 5 ms to 20 ms (to omit initial transient effects). The power spectral density 
of u(xj,t) at each j is obtained using signal processing tool (sptool) in Matlab, to illustrate the 
energy distribution across frequencies in log-scale. 

Figure (jj) top frame shows the time series of BM displacement at x = 1.93 cm, the bottom 
frame is the log-log plot of FFT power spectral density vs. frequency. We see the major peak at 
3.5 kHz, as x = 1.93 cm is the (so called characteristic) location for the peak of a single 3.5 kHz 
tone. In addition, we see two small side peaks at 3 kHz (= 2f\ — / 2 ), and 4 kHz (/ 2 ). In figure 
at x — 1.85 cm, the characteristic location for / 2 , the / 2 peak is more pronounced, however 
the fx peak is still the highest. Such effect of lower frequency tone fx to a higher frequency tone 
/ 2 is called upward masking in hearing. In addition, there are two small side peaks at 3.0 kHz 
( = 2 fx — / 2 ) and 4.5 kHz (= 2/ 2 — fx). In figure (j3J), at x — 2.03, the characteristic location 
for 3 kHz, a dominant single peak due to the generated combination tone 2fx — / 2 is observed. 
In contrast, the 2/ 2 — fx tone (= 4.5 kHz) is weaker and dominated by fx and / 2 even at its 
characteristic location x = 1.69 cm, see figure (J3J). The above findings on combination tones 
are consistent with the experiments on cochlea and the analytic structures of multi-tone 
solutions in the previous section. 
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5 Conclusions 



The nonlinear nonlocal cochlear models of the transmission line type are well-posed globally 
in time and admit exact multi-frequency solutions in the weakly cubic nonlinear regime. For 
finitely many tonal input at distinct frequencies, the exact solutions contain all integral linear 
combinations of input frequencies. For two tone input with frequencies fx and f 2 at high 
enough intensities, we observed numerically the combination tones 2fi — f 2 and 2f 2 — fi in 
model output, in agreement with existing experimental observations [T7j and the structure of 
analytical solutions. 
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BIVI time response at location x — 1 .93 cm 




"Time in ms 

FFT power spectral density at BM location x=1 .93 cm 



20 



-40 



-60 



-80 



I 

Input signals: 80 dB 3.5 kHz 


i 


i 


85 dB 4.0 kHz 










*• 3.5 kHz 




3 kHz ^ 

















10 1 10 2 10 3 10 4 10 5 

Frequency (log-scale) 



Figure 1: Time series of BM displacement at x = 1.93 cm (top frame), and its FFT power 
spectral density vs. frequency (bottom frame). 
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Figure 2: Time series of BM displacement at x — 1.85 cm (top frame), and its FFT power 
spectral density vs. frequency (bottom frame). 
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BIVI ti me response at location x= 2.Q3 cm 
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FFT power spectral density at BM location x= 2.03 cm 



Input signals: 80 dB 3.5 kHz 
85 dB 4.0 kHz 



10 



10' 



10° 

Frequency (log-scale) 



3 kHz 



10' 



10" 



Figure 3: Time series of BM displacement at x — 2.03 cm (top frame), and its FFT power 
spectral density vs. frequency (bottom frame). 
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BIVI time response at location x - 1 .69 cm 
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Figure 4: Time series of BM displacement at x = 1.69 cm (top frame), and its FFT power 
spectral density vs. frequency (bottom frame). 
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